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We present a phase diagram for a dilute two-dimensional Bose-gas on a lattice. For one sort 
of boson we consider a realistic case of the van der Waals interaction between particles with a 
strong hard-core repulsion U and a van der Waals attractive tail V. For V < 2t, t being a hopping 
amplitude, the phase diagram of the system contains regions of the usual one-particle Bose-Einstein 
condensation (BEC). However for V > 2t we have total phase separation on a Mott-Hubbard Bose 
solid and a dilute Bose gas. For two sorts of structureless bosons described by the two band Hubbard 
model an s-wave pairing of the two bosons of different sort (6162) 7^ is possible. The results we 
obtained should be important for different Bose systems, including submonolayers of 4 He, excitons 
in semiconductors, Schwinger bosons in magnetic systems and holons in HTSC. In the HTSC case 
a possibility of two-holon pairing in the slave-bosons theories of superconductivity can restore a 
required charge 2e of a Cooper pair. 
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I. INTRODUCTION 



It is well known that, in contrast with two-particle Cooper pairing in Fermi systems, the essence of a superfluidity 

in Bose systems is one particle Bose-Einstein condensation (BEC). This asymmetry between FermLand Bose ( two- 

~~ . particle versus one-particle condensation) was challenged in a pioneering paper by Valatin and ButlerEJ. They proposed 

i__i' a BCS-like variational function for the description of an attractive Bose gas. The most difficult problem with the 
validity of their description is connected witlua tendency toward phase separation which arises in attractive Bose 
systems. Later on Nozieres and Saint- JamesEl conjectured that in a Bose system with a short-ranged, hard core 
repulsion and a van der Waals attractive tail, in principle, it is possible to create a two-particle bosonic bound state 
and to escape collapse. Unfortunately their calculations in a three-dimensional system showed that at least for one 
sort of structureless boson either standard one-particle BEC is more energetically beneficial, or that phase separation. 

(^. . takes place earlier than the two-particle condensation. Note that the same result was obtained earlier by Jordanskicl 

<^> for the case of weak van der Waals attraction. n 

The important development of the ideas of Nozieres and Saint James belongs to Rice and Wangu. The authors 
claimed that in two-dimensional (where already an infinitely small attraction leads to the bound state in a symmetrical 
potential well) it is possible to realize a two-particle boson pairing. Moreover, this two-particle pairing results, for 
small momenta p£o < 1 in a linear, soundlike, dispersion law of quasiparticles at T — in an analogy with a standard 
one-particle Bose-condensation. 

f^j . To escape a collapse in a 2D attractive Bose gas, the authors of Ref. introduced in their model a Hartree-Fock 

shift of the chemical potential [Ib ~ Un, connected with the short-range repulsion U. This shift in the case of 
U > 2V, where V is the magnitude of the van der Waals tail, leads to a positive compressibility in the system 
k- 1 = dus/dn = U - 2V > 0. 

The main goal of a present paper is to construct a phase diagram of a 2D dilute Bose gas with the van der Waals 
interaction between particles, by taking into account on equal grounds the full contribution of a hard-core repulsion 
U and a van der Waals tail V. Throughout the paper we will consider the lattice model, and will base our results on 
the exact solution of the two-particle T-matrix problem presented in Refs. @ and a We will study the possibility of 
different two-boson pairings, as well as the possibility of a total phase separation in the system. We will also consider 
the two sorts of structureless bosons described by the two-band Hubbard model. In the case of attraction between 
bosons of two different sorts, we will find a possibility of an s-wave two-boson pairing (6162) 7^ 0. 

II. THEORETICAL MODEL 

The model under consideration is described by the following Hamiltonian on the 2D square lattice: 



ij i ij 

where n, = 6j6^ is a 2D boson density. We will work in the limit of strong hard-core repulsion U 3> {V; t}, and restrict 
ourselves mostly to a low-density limit ngd 2 <C 1, d being the interatomic distance. Further on we put d = 1. 

Note that in the case of V = model (H) is just the Bose-Hubbard model, extensively studied in the literature for 
the case of 2D 4 He submonolayers, as well as for the flux lattices and Josephson arrays in the type-II superconductors 
(seeRefs. g— fu]). 

We would like to emphasize that model ([jj) is, to some extent, a Bose analog of the fermionic t — J model with 
released constraint, considered by Kagan and RiceQ. After the Fourier transformation from Eq. (ffl) we obtain: 

P kik 2 q 

- E ^(q) & kAiW>kA 2 - q , (2) 

kik 2 q 

where e p = — 2t(cosp x + cosp y ) is a bosonic spectrum on the square lattice and V(q) = V(cos(q x ) + cos(q y )) is a 
Fourier transform of the Van der Waals tail. As a result, a total interaction in the momentum space is given by the 
formula: 

Veff(q) = ^-V(q). (3) 



III. THE T-MATRIX PROBLEM. 

An instability toward a two-particle boson pairing manifests itself in the appearance of a pole at a temperature 
T = T c in the solution of the Bethe-Salpeter equation for the two-particle vertex T for zero total momentum of the 
two bosons (p, — p). To proceed to the solution of this equation, we must solve at first the T- matrix problem for 
the two particles in vacuum. Here we can use the results of Ref. g, because the solution of the two-particle problem 
does not depend upon statistics of colliding particles. For the T- matrix problem it is convenient to expand V e ff(q) 
in Eq.(S) in series with the eigenfunctions of the irreducible representation of the lattice symmetry group -D4. This 
yields 

V e ff (extended s-wave) = ^ — y( cos Px + cosp y )(cosp' x + cosp' y ), 

V e ff (p-wave) = -^-(smp x smp' x + s'mp y smp' y ), (4) 

V eff (d x 2_ y 2) = -¥-(cosp x - cosp y )(cosp' x - cosp'y), 

where q = p — p' is a transfered momentum. Note that, for spinless bosons, which we formally consider in Eq. (1), 
the total spin of the Bose pair is zero. Hence only s- and a d-wave pairings are allowed by the symmetry of the pair \P 
function. A p-wave pairing is allowed only for an odd total spin of the two bosons. Nevertheless we will conserve the 
results for the p-wave pairing in our paper because the generalization of Eq.(l) for the case of bosons with internal 
degrees of freedom is straightforward. The T-matrix problems for p- and d-wave channels are very simple. Solutions 
of these problems for the two particles with a total momentum zero and a total energy E yield 

T d -JE) = --, \ r , (5) 

where 

= f f dp x dp y \c/) d , p \ 2 

d ' p J J (2ir) 2 E + 4t(cosp x + cosp y ) , {> 



and the functions <j) d — cosp x — cosp y and (j) p = smp x + i smp y correspond, respectively, to the d- and p-wave channels. 
Let us find the thresholds for the bound states in the d- and the p-wave channels. The appearance of a bound state 
means that E = —W — \E\, where W = Bt is a bandwidth. For the threshold E = 0. An exact solution of Eqs. (H) 
and (ph, which involves the calculation of elliptic integrals of first and second order, yields 
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Note that a threshold for a p-wave pairing is lower. Now let us proceed to an s-wave channel. Here an ordinary 
s-wave pairing is suppressed by large hard-core repulsion U, however an extended s-wave pairing with a symmetry of 
the order parameter A s = Ao(cospj; + cospj,) is allowed. In real space this pairing corresponds to the particles on 
the neighboring sites. Moreover the pair \P function is zero in the region of a hard core (r < ro) and is centered (has 
a maximum) in the region of a van der Waals attraction for r ~ r\ (see Fig. 1). On the lattice ro ~ d/2 and r% ~ d. 
One can see that the ^> function has a region of zero value zero. But it has no nodes because it does_npt change its 
sign for all values of r. The rigorous calculation of the threshold for an extended s-wave pairing givesEKa 

I) - 1- < 8 » 

^ b / s-wave 
Moreover for V > V c an energy of the bound state has the form 

w-fl-^tav^b}- (9) 

Of course, in a strong coupling case V 3> W 

\E b \ « V. 
The T matrix in an s-wave channel for small and intermediate values of V is given by 

r -(' £ i> - w n ln g v - < 10 ' 

The most important is that a strong Hubbard repulsion U acts only as an excluded volume and effectively drops out 
from Eq. (^Tj) at low energies. For V <§C At the T matrix 



r.(|£|) 



TtW 



corresponds to repulsion and coincides with the T matrix for the 2D Bose-Hubbard model at low density. Of course, 
at very high energies the T matrix contains an additional pole E rs U which corresponds to an antibound state in a 
total analogy with the fcrmionic Hubbard model. For V — At: T S (\E\) = and there is no interaction at all. Finally, 
for V > At, T S (\E\) < corresponds to an attraction and reflects the appearance of the bound state. 

IV. BETHE-SALPETER EQUATION FOR AN S- WAVE PAIRING OF THE TWO BOSONS 

Let us consider at first the most interesting case of V > At and find the critical temperature for jail extended s-wave 
pairing of the two bosons. The solution of the Bethe-Salpeter equation for bosonic systems readalJ 

T 

r s = r— • (ii) 

" dp x dpy coth £ ^ T P 



1 + T t 
For a low density of bosons tib <S 1 one has 



(2tt)2 2(e p -/i) 






-p 



where m = l/2t is a boson mass. Accordingly /i = — At + jj, and £ p = e p — u = p 2 /2m — jl. 

The most substantial difference of Eq. ( |ll| ) from an analogous fermionic equation is the replacement of tanh £/2T 
by coth£/2T in its kernel. Moreover, as shown in Ref. [lj for the 2D Fermi gas jl — Ep~ \E\\j1. So, in a weak-coupling 
case, when ep ^S> \Ej,\, the chemical potential jl « ep > is positive. In contrast to this we shall see below that a 
bosonic chemical potential jj, is always negative even in the weak-coupling case, when a binding energy is much smaller 
than a degeneracy temperature: 

|£ 6 |«T = ^. 
m 

Another very important point is that the T matrix, which enters into the Bethe-Salpeter equation, must be calculated 
for a total energy E = 2/1 of colliding bosons. The chemical potential jl can be determined from the requirement of 
the number of particle conservation. This requirement yields 

f f d 2 p 1 

nB = iiWe^^}TI- (12) 



From Eq.(|12j) for the temperatures \Eb\ < T < T < W we obtain 

A = -Texp(-^ <0. (13) 

Note that a standard Hartree-Fock shift nil drops out from the expression for g p = e p — u both in the Bethe-Salpeter 
equation (O) and in the equation for the number of particle conservation (12). Now we are ready to solve the 



Bethe-Salpeter equation (|ll|). The critical temperature T c corresponds to the pole in Eq. (Jill 



mTJ2a) 

1 + 2 y J I = 0, (14) 



where 
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An analysis of Eq. ( 19) shows that the main contribution to the integral comes from the lower limit of integration. 
Hence providing \p, 



/T c «lwe have 

W/T c 

<k 



,.. i dy _ 2T C 



!! ,y '2T, 



As a result Eq. ( pT| ) can be represented in the following form: 



\p\ mT a {2p) 

It is useful now to represent T s (2jj) in terms of the binding energy Ef,. Utilizing Eqs. d|) and (00) we can write 



(17) 



Ts (2~a) = -J^L = *!_ (18) 

lnM m l n M 



r^r " ,m iB 



It is important to mention here that ft < 0, and hence the T matrix in Eq. (18) does not contain an imaginary 
part. In the fermionic case ft = Ef > 0, and the T matrix contains an imaginary part corresponding to the resonant 
scattering. As a result, from Eq. (J18|) we obtain: 



41 °f = !f- < w > 



Assuming that \E b \ <C T c <C To , we obtain 



and 



ft{T c ) = -T c cxp ( ~y 



IAI . f T Q 



Later on we will justify this assumption. 
As a result from Eq. ([H]) we will obtain 



t= **»}-¥„ 



T * = \ ( ^ m y (20) 



Recall that in the case of the fermionic pairing in two dimensionsa critical temperature readsc£l 

T c = ^2e F \E b \. 
Let us analyze expression (|20|). As we already know 



where 



Then a condition \E b \ <Tp means 



I./IV ■ Mr.-. M ><i-~ 



A=<^>. (21, 



A« T -V«1. (22) 



ln-^- 
111 To 



Hence \nT /\E b \ = 1/A - In W/T w 1/A, and 



U1 4 1U |B b | m 4 A 



which is in an agreement with Ref. ^. Note that T c from Eq. (|23|) satisfies the conditions \Eb\ <C T c <C To, so an 
assumption used for the derivation of T c is justified. 

For T < T c the spectrum of the quasiparticles acquires a gap: 



E P 



^/( £p + |A|) 2 -A2. (24) 



Note that at low densities of bosons a gap A becomes isotropic in the principal approximation. 

The gap A together with the chemical potential ft must be defined self-consistently from the two coupled equations 



,"_ 



X f coth Jz 2 - -^r 
1 = 5 / dz / (25) 

6 4^T 



nB = I / d t f J—- 1 > ( 26 ) 

IAI exp|^ j-1 

where £ = £ + \ji\ and 2 = £/2T. 

Of course, the solution of the system of equations Eqs. ( p5[ ) and ( |2q ) exists only if |/2| > A, or, in other words, only 
if Eg — \p\ 2 — A 2 > 0. The exact solution of these equations yields for zero temperature in an agreement withQ 

|A(T = 0)|=A=J^1. (27) 

This result is very important. It justifies our scenario, leading to a linear, soundlike spectrum of the quasiparticles 
for a small momenta p. Indeed 



E p = ^sl + 2e p \tt = ^el+e p \E b \. (28) 

From Eq. ( p8|) for the case p£o <C 1, where £0 = l/y/2m\Ei,\ is the coherence length of the boson pair, we immediately 
obtain a linear dispersion law: 

E p = cp. (29) 

In ( |29| ) c 2 = \Eb\/2m is a sound velocity squared. This means that an inverse compressibility of the system k^ 1 = c? 
is positive. This fact proves the stability of a superfluid paired state and excludes the possibility of the collapse of 
the pairs in the system. Note also that close to T c one has 



A(T)«A(0)^p^, (30) 

which is similar to the BCS theory. We would like to mention that bosonic pairs in the limit \Ei,\ <C T are extended 
in full analogy with the BCS theory. That is the coherence length in this limit, 

Co » ^= » 1, (31) 

is larger than the mean distance between the bosons. The Bose pairs are strongly overlapping in this limit. The pairing 
takes place in the momentum space in an analogy with the Cooper pairing in the BCS picture of superconductivity. 
In the opposite limit \E^\ ^> T , the situation closely resembles the bipolaronic limit for thej-iermionic systemsEao. 
That is, the creation of the bosonic bound pairs is associated with the crossover temperaturefij 

The Bose condensation of the pairs occurs at lower temperature 

lnln(l/n) 

Note that this temperature is obtained from the condition fib(T c ) = — T c cxp (— To/T c ) + foT c = 0, where fa = 
1/ ln(l/n) is a repulsive interaction between the local pairs. So the superfluid transition takes place only for a residual 
repulsion, between the pairs. Also note that in a dilute Bose gas the Berezinski-Kosterlitz-Thouless contribution of 
vorticeala is important only very close to T c , so the mean field expression ( p3|) gives a very good estimate for the exact 
BKT critical temperature: 

{T c - T BK t)/T c ~ l/lnln(l/n) « 1. 

In the case of the local pairs the coherence length is small: 

Co « 4=- ( 34 ) 

The pairs are compact, and the pairing takes place in the real space. 

G 



V. POSSIBILITY OF P-WAVE AND D-WAVE PAIRING OF THE TWO BOSONS 

Now let us analyze the solution of the Bethe-Salpeter equation for p- and d-wave boson pairings. Here the critical 
temperatures should be found from the conditions 

1 + T pid (2ji)ir,d = 0, (35) 



where 



dp x dp y coth -^r- | ^ 
(2tt)2 2( £ - M ) 








In a low-density limit the (f> functions can be approximated by the following expressions: 

(j> P ~Px + ip y = pe l<t ', 



h ~ ~(p£ - p\) = -P 2 cos 20. 



Hence after an angular integration we obtain: 

to f coth ^ 

i p = ^z P d P — w^ p ' 

7 r coth £ (37) 

where again £ = p 2 /2m + \j2\. 

Additional factors p 2 and p 4 in the integral expressions for I p and 1^ reflect a well-known fact, that for slow 2D 
particles in vacuum an s-wave harmonics of the scattering amplitude behaves as /o ~ In 1/p 2 , whereas for a magnetic 
number m ^ 0, the scattering amplitude vanishes for p goes to zero as f m ~ p 2m . The additional factor p 4 leads to 
the absence of an infra-red singularity for e — > in 1^: 

f dee 1 
h ~ / — =- ~ e -► 0. (38) 



For the p-wave channel the infra-red singularity becomes logarithmically weak: 

J p ~y^~lne. (39) 

This means that the Bcthc-Salpeter equation has no solutions in p- and d-wave channels for |V|/i < 1. 

Hence the boson pairing with a large coherence length £o > 1/ ' ^fn is absent in a p-wave channel as well as in a 
d-wave channel. Here only the limit of the local pairs is possible. For p- and d- wave channels local pairs are created 
at the crossover temperatures given by Eq. (B2). In this case the binding energies for p- and d-wave channels read, 
correspondingly, 

\K\*W {V - Vcp) \ , (40) 

\Et\ « W iV ~ v Vcd) . (41) 

Note that V cp and V c d in Eqs. (40) and (|4l| ) are the thresholds for the p- and d-wave pairings given by Eq. (0). We 
would also like to mention that for a fixed V: \E%\ > \E*\. Providing that the interaction between the local pairs is 



repulsive, the temperature of the Bose condensation of the local pairs in these two cases is again given by Eq. (pj 
Summarizing this case we see that, while the crossover temperatures T* are different for the s-, p- and d-wave local 
pairs, their critical temperatures coincide with each other in a principal order. 



VI. TOTAL PHASE SEPARATION 

As we discussed in the Sec. 1, the real collapse is prohibited in our system by large Hubbard repulsion U. However 
the phase separation on the two large clusters is allowed. The first cluster corresponds to the Mott-Hubbard Bose 
solid. In this cluster n& — ► 1, that is, each site is practically occupied by one boson. Such a cluster is localized due to 
Mott-Hubbard considerations. It has no kinetic energy. However it has a potential energy of the order of — 2V for one 
particle. A second cluster has a very small boson density n — > 0. In this cluster for V < At the energy per particle is 
— W/2 + fan, where /o = m/47rln 1/n is a two-particle T matrix in the absence of a bound state. Rigorously speaking, 
at a given bosonic density n the phase separation results in the formation of the two clusters with the densities n\ > n 
and ri2 < n, where n\ is close to or identically equal to 1. The phase separation for V < At takes place if the energy 
per particle in the cluster with the density n\ becomes smaller than the energy per particle in the cluster with the 
density ri2 

- 2j]V < -W/2, (42) 

where r\ is an unknown numerical coefficient of the order of 1. Note that in the fcrmionic t — J model, considered in 
Ref. [| the Mott-Hubbard cluster with n — 1 has an antiferromagnetic order. Hence instead of T]V in ([42|) one should 
write 1.18J - the energy per spin of the 2D antiferromagnetic. As a result, in a fermionic case J ps — 3.8i. In our 
system V ps ~ 2t, due to the absence of kinetic energy and zero point energy in the case of structureless bosons. In the 
same time for n — ► 1 the phase separation between the Bose solid and the one-particle BEC takes place. According 
to Ref. |l^ for Hb — > 1, the phase separation takes place already for small values of |V|/t. r— ■ 

In principle, another scenario of the phase separation connected with the creation of quartetsEj is also possible in 
our system. It requires an evaluation of the four-particle vertex which is impossible to do analytically. However we 
think that our scenario of total phase separation-takes place for smaller values of V/t than the quartet formation. 
This is in agreement with numerical calculationsta for the 2D fermionic t — J model. 

VII. PHASE DIAGRAM OF THE SYSTEM 

In this section we will complete the phase diagram of the system. At first, note that for V < V ps (when the T 
matrix for an s-wave channel is repulsive) we have at low density a standard Bose gas with a hard-core repulsion. It 
will be unstable toward a standard one-particle BEC at a critical temperature given again by Eq. (|3^). For V > 2t a 
total phase separation on two large clusters takes place in our system. One of these clusters contains a Mott-Hubbard 
Bose solid, another one contains a Bose gas with one particle condensation. 

For large densities n — n c < 1 (n c = 1 in Ref. [lOJ for structureless bosons) the system will undergo a transition 
to the Mott-Hubbard Bose solid. As a result, on a qualitative level the phase diagram for our system has the form, 
presented in Fig. 2. Note that our phase diagram should be important for the understanding of the physics of the gas 
of kinks and steps on a solid interface of 4 He. It will be interesting in this context to obtain phase separation in a 
system of kinks on a quantum atomically rough surface of 4 He. Another possible application of our results can be 
connected with the study of biexcitonic pairing in semiconductorsEj. In this context we must mention an interesting 
recent paper of Lozovik et al .123. 

VIII. TWO-BAND HUBBARD MODEL FOR THE TWO SORTS OF BOSONS 

Let us consider the two-band Hubbard model for the two sorts of structureless bosons. The Hamiltonian of the 
system has the form 



H = -t a Y^ a l a ~ *«• 2 ^ b i 



y y 

+ ~?rJ2 n ™ + -!fJ2 n rt- -^J2n ta n lb , (43) 

i i i 

where t a and tb and n a and rib ar e, respectively, the hopping matrix elements and densities for bosons of sorts a and 
b. For simplicity we will consider the case t a = tb, corresponding to the equal masses m a = l/(2t a ) = nib. We also 
assume that the bottoms of the bands coincide. In the Hamiltonian U aa and Ubb are Hubbard onsite repulsions for 
bosons of sorts a and b. Finally U a b is an onsite attraction between bosons of two different sorts. 



Let us consider the low-density limit, when both n„ < 1 and n b -c 1. In this limit we must replace the Hubbard 
interaction U ab by the corresponding T matrix. The relevant expression for the T matrix T ab is given by 

T *b(E) = : TT r U l i (44) 



1 - U ab J 



d 2 i> 



4-zr 



p 2 /m+\E\ 



where E is given again by E = E + W . The T matrix has the pole for the energy 

\E\ = \E b \=W^v{--^-\. (45) 

In the extremely strong coupling case U a b > W the pole corresponds to the energy |J5| = \E b \ = U ab . The pole in the 
T matrix reflects the appearance of the bound state of the two bosons of different sorts. 

Now we can solve the two-particle problem in the presence of the bosonic background. A simple analysis shows 
that only local bosonic pairs (bipolarons) are possible in our case. They are formed at a high crossover temperature 
T* = \E b \/ln(W/T Q ), where T = min{To ,T 6 }, and T a and T 6 are degeneracy temperatures for bosons of the two 
sorts. Correspondingly the pairs are bose condensed at lower temperature T c ab = T /ln(l/A), where A = mU a b/{^)- 
Our results are valid is the case \E b \ > {Tq ,Tg}. In the opposite case of higher densities, when at least one of the 
temperatures Tq or Tq is larger then \E b \, we have at first a standard one particle condensation for bosons with higher 
density. As a result the two-particle pairing between bosons of different sort can take place only as a second superfiuid 
transition inside the superfiuid phase with one particle BEC. Let us now analyze the stability of our system with 
respect to collapse. For simplicity we consider an extremely strong coupling case {U aa , U bb , U a b} > W. In this case 
the local pairs of two bosons of different sorts have onsite character (aj&i) ^ 0. To escape collapse in this case we must 
satisfy the following stability criteria: namely, we must prevent the collapse of (ab) pairs on one site (a creation of the 
quartets (aabb)). This requires the inequality U aa + U bb — ^U ab > 0. These results can be important for slave-boson 
and Schwinger-boson theories. 

IX. SLAVE-BOSON FORMULATION OF THE T- J MODEL. 
APPLICATION TO HTSC. 

We conclude our paper with a brief discussion of the problems which arise in a slave-boson formulation of the t -*-J 
model. In this formulation close to half- filling (n — > I) an electron is represented as a product of holon and spinoncJ: 






= flbi 



A superconductive d-wave gap 

\i = (c^Cji - c iiCn ) 
is a direct product 

A d = A sp A{ 
of a spinon <i-wave gap 

A sp = (fitfjl -/u/it) 
and a holon s-wave gap 

A fc = <6 i 6 i ). 

Then a natural question arises whether (bi) ^ and, accordingly, A/j = (bi)(bj), or (b t ) = but (6.^^=^ 0. In other 



words, whether a one particle or two-particle condensations of the holons takes place in our systemol 

This problem is a very difficult one and, surely deserves a very careful analysis. Our preliminary considerations show, 
however, that the mare beneficial conditions for the two-particle condensation may arise in the SU(2) formulation 
of the t — J modelESED. In the standard U(l) formulation of the modeled an effective potential of the two-holon 
interaction on neighboring sites appearing after the Hubbard-Stratonovich transformation has a form 



T-l)XW*4* (46) 

ij 

and thus corresponds to the repulsion for t > J. This observation excludes the possibility of the two-holon pairing in 
the U(l) formulation of the t — J model. 

In the SU(2) case it will be desirable to derive conditions when (b\) = (62) = but (6162) 7^ 0. For such a 
nondiagonal pairing, as already discussed above it is easier to satisfy the stability criteriaa. Note also that the same 
situation with two sorts of bosons and a possible attraction between them can be realized for 2D magn etic systems. 
We obtain the corresponding bosonic Hamiltonian here after a Schwinger transformation of spinsE2rc3 in extended 
Hcisenberg models. Work along these lines is now in progress. 

X. CONCLUSIONS 

In conclusion we analyzed the possibility of the formation of boson pairs with s-wave symmetry and an appearance 
of total phase separation in a 2D Bose gas. In addition we considered the case of boson pairs with symmetry of p- 
and d-wave type. We also constructed the qualitative phase diagram for the 2D Bose gas with the van der Waals 
interaction between the particles, which, besides a standard one particle BEC, contains regions of the Mott-Hubbard 
Bose solid and a total phase separation. We also consider the situation for two sorts of bosons described by the 
two-band Hubbard model, and found the conditions for the two-particle pairing between bosons of different sorts. 
We discuss the applicability of our results for the different physical systems ranging from submonolayers of 4 He and 
excitons in the semiconductors till Schwinger bosons in magnetic systems and holons in the slave-boson theories of 
HTSC. 
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FIG. 1. The ^-function of an extended s-wave pairing, ro is a radius of a hard core repulsion, r 1 is a radius of a Van der 
Waals attraction. On the lattice ro ~ d/2 and n ~ d 
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FIG. 2. Qualitative phase diagram of the 2D Bose gas with the Van der Waals interaction on the square lattice. 
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